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CONTACT CENTER SCHEDULING USING INTEGER 
PROGRAMMING 

FIELD OF THE INVENTION 

The invention pertains to the field of computer systems and computer 
implemented methods. More specifically, the invention pertains to methods for 
scheduling agents with multiple skills optimally at contact centers to satisfy time varying 
agent staffing requirements for multiple contact types. 



BACKGROUND OF THE INVENTION 



The need for scheduling agents optimally in contact centers arises in a broad 
5 range of industry areas. Examples include call centers providing inbound or outbound 
services, telephone operator centers, catalog sales centers, hospitals, police stations, 
airports, maintenance and installation operations, manufacturing, etc. Agents may be call 
center operators, technicians, workers, servers, nurses, etc. Contact centers may include 
any service center providing services to its customers through their agents. Many 

10 manufacturing operations may also fit into a contact center paradigm with similar needs 
and challenges in terms of worker or crew scheduling. Scheduling agents, may be done a 
day at a time by specifying daily shift start time, shift length, time for the relief and lunch 
breaks, or to cover an entire week specifying, in addition to specifying daily shift 
schedules, work and non-work days. Contact centers use one or more interaction media 

1 5 to provide the desired services to their customers. Interaction media (also referred to as 
"contact types") may include telephone calls, web chat^ email, fax, etc. It may also 
include contacts using different languages such as English, Spanish, etc., and/or different 
services such as sales, technical help, etc. Customers choose their time and media of 
contact to interact with the contact center agents. The contact type used by a customer 

20 may require an agent with a special skill (e.g. language skills, writing proficiency). 

Moreover, agents may possess only some of the skills needed to serve customers using 
different contact types. Agentis with the same skill set form a skill group (also referred to 
as "agent skill group''). 

25 One of the critical tasks in contact center management is workforce scheduling. 

In such systems, workload changes during the course of an operating day, and fi-om one 
day to another during a week. While understaflBng will lower labor costs, it will result in 
poor service levels and longer waiting times for customers. Overstaflfing, on the other 
hand, will improve service levels and lower waiting times, but will result in underutilized 

30 agents and higher labor costs. Hence, it is important for contact center managers to 

schedule their workforce in an optimal manner to minimize labor costs while meeting the 
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projected agent and skills requirements. The present invention provides a method for 
formulating the workforce scheduling environment and requirements at a plurality of 
contact centers as a Mixed Integer Linear Programming (MILP) model, and a solution 
algorithm for developing optimal schedules for agents. 

5 

Traditionally, the task of developing work schedules to meet the forecasted agent 
and skill requirements is carried out manually. Manual scheduling when there is time 
varying workload and agent requirements, various daily shift and weekly tour scheduling 
rules, and multiple agent groups with different sets of skills and full-time/part-time work 
1 0 requirements is laborious and inefficient. Manpower scheduling in a typical contact 
center environment may involve an astronomical number of potential schedules, making 
manual evaluation of even a very small fraction of them to find an efficient schedule 
impractical. 

15 In the past, a number attempts were made to develop more effective methods for 

workforce scheduling. Even though terms such as "optimal schedule", "optimized 
schedule" or "local optimal" were used in their descriptions, none of the prior art agent 
scheduling methods states or verifies the necessary and sufficient conditions for 
optimality in scheduling agent. Like the manual scheduling method, the prior art agent 

20 scheduling methods use a set of rules of thumb to develop agent schedules. None of 
them presents any evidence of, or guarantee the optimality of the schedules that may be 
obtained using the rules of their methods. 

In the field of mathematical optimization, there are proven necessary and 
25 sufficient conditions for a solution to an optimization problem to satisfy in order to be the 
optimal (absolutely the best) solution (Taha, 1987, pages 743-779). The method of the 
present invention first develops an MELP model for which these necessary and sufficient 
conditions are well stated, and then uses an algorithm that would generate an optimal 
solution satisfying these conditions. 

30 
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U.S. Pat. No: 5,91 1,134 issued on June 8, 1999 to Castonguay et al. describes a 
method for developing weekly tour schedules in call centers for agents with a single skill 
type. This method describes tour and shift construction rules to select tour and daily shift 
templates from the available ones, one day at a time, using measures such as "coverage", 
5 and break scheduling rules to schedule daily breaks (U.S. Pat. No. 5,911, 134, col. 18, 
lines 26-46, and FIG. 1 1). The method also includes post-processing steps to ehminate 
the redundant tours, level break times, and take employee preferences into consideration. 
No evidence to suggest that the schedules obtained using these rules will satisfy the 
necessary and sufficient conditions for optimality in a workforce scheduling problem is 
1 0 provided or known. As in the case of manual scheduling, the existence of an 
astronomical number of potential schedules, even in small call center scheduling 
environments, makes it very unlikely to locate the optimal schedule with the tour and 
shift construction rules used by this method. 

15 U.S. Pat. No. 6,044,355 issued on March 28, 2000 to Crockett et al. describes a 

simulation method for developing weekly tour schedules in a contact center environment 
involving multiple agent groups with a plurality of skill sets, and a plurality of contact 
types requiring different agent skills. The method uses a scheduler and an Automatic 
Call Distributor (ACD) simulator. However, no working details or description of the 

20 scheduler, and how it develops a schedule is included in the patent document. It was 
disclosed that the method uses a "scheduler program" (U.S. Pat. No. 6,044,355, col. 6, 
lines 18-36). No evidence to suggest that the schedules obtained using the scheduler 
program considered will satisfy the necessary and sufficient conditions for optimality in a 
workforce scheduling problem is provided or known. 

25 

U.S. Pat. No, 6,278,978 issued on August 21, 2001 to Andre et al. describes a 
method that may be used to post-process an initial schedule available from another source 
outside of the method by applying a rule-based interchange procedure. The method is 
provided for scheduling agents with a single skill type. The method begins by acquiring 
30 an agent schedule developed by another method. The method then unschedules an agent 
from a shift in the schedule, evaluates the value of a "score" fiinction for all other 
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schedules that may assign this agent to a different shift, and selects a shift using the score 
function to improve the schedule. The method also considers rescheduling of breaks, 
again using the score function. Interchanging two or more activities to improve the 
quality of an initial solution to a problem is a well-known strategy in the field of 
5 Optimization and the related research literature. It is also well known in the field of 
Optimization that interchanging (rescheduling) one or more agents and shifts 
simultaneously doesn't guarantee the optimality of the solution developed and the 
satisfaction of the optimality conditions. 

A primary object of the present invention is to overcome these limitations of the 
prior art agent scheduling methods. The present invention uses Mixed Integer Linear 
Programming approach (IVDLP) (Wolsey, 1998) to develop a mathematical model of a 
workforce scheduling environment. A number of researchers including Danzig (1954) 
proposed MILP models for workforce scheduling. Difficulties with the use of the prior art 
MILP models for workforce scheduling are well documented ( Nanda and Browne, 1992, 
page 206, Holloran and Bym, 1986, page 13). The MILP model proposed by Dantzig 
(1954), for example, enumerates all possible combinations of shift and tour parameters 
resulting in tens of trillions of decision variables, making it very inefficient and 
impossible to solve even in small contact center environments. To overcome these 
limitations of the prior art MBLP agent scheduling models, the present invention 
formulates daily break, and days-off scheduling implicitly (i.e. not explicitly enumerating 
all possible combinations). Thus, the method of the present invention formulates a 
significantly smaller but equally powerful MILP model, and uses various extensions of it. 

25 A ftirther limitation of the prior art MILP models is that they only consider agent 

scheduling environments involving a single agent skill group. Recent technological 
developments in ACD's used in call centers to queue and assign calls to agents, and 
incorporation of other types of contact media such as email and fax, made contact center 
managers to realize that agents with different types of skills are needed to be scheduled to 

30 handle these contact types (e.g. a Spanish speaking agent to answer a caller speaking in 
Spanish, or agents with good writing skills to respond to emails). Thus, the agent 
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scheduling task became more complex since the agent scheduling method used should 
take the time-varying demand for different agent skill types, and available agent groups 
with different skill sets (e.g. Spanish & English speaking sales, Chinese & English 
writing for technical help emails, and English only agent groups) into consideration in 
5 scheduling agents. It is another object of the present invention to disclose a method for 
formulating MILP models for scheduling environments involving a plurality of agent 
skill groups, and a plurality of contact types with time-varying workload for specific 
agent skills. 

10 The present invention also discloses a solution algorithm that solves a number of 

sub-problems (or nodes) with the use of the standard Branch and Cut (B&C) algorithm 
for MILP problems. The solution method supplements the B&C algorithm by a 
Rounding Algorithm (RA algorithm) to locate the optimal solution of the MELP model of 
the invention. The optimality condition (both necessary and sufficient) for an MILP 

1 5 model is well documented in the field of Optimization; An all-integer solution that 
satisfies all of the constraints of an MILP with a minimization (maximization) type 
objective function, and has an objective value that is better than the best lower bound 
(upper bound) for any of the sub-problems in the B&C algorithm is a global (i.e. 
absolute) optimal solution to that MELP model (Wolsey, 1998). Once the optimality 

20 condition is satisfied by an all-integer solution found during the execution of this solution 
algorithm (either in the B&C algorithm or in the RA algorithm), an optimal agent 
schedule is reported using the information in the optimal solution of the MILP model 
found. 



25 It is still another object of the present invention to provide a computer 

implemented optimal workforce scheduling method for contact centers that schedule their 
workforce to meet varying workload during the course of a day, and fi'om one day to 
another during a week. 

30 
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SUMMARY OF THE INVENTION 



The method of the present invention first formulates a skills based or non-skills 
based agent scheduling environment involving a plurality of agent skill and work groups, 
5 and plurality of contact types as an MILP model including various contact center 
constraints and requirements. The formulation step of the present invention involves, 
first, the definitions of various decision variables, and parameters, and the formulation of 
an objective fiinction. The objective Sanction of the MELP model of the present invention 
determines a measure of effectiveness (or merit) for each solution. Examples of the 

10 objective Sanction in the method include total schedule cost, total time scheduled, total 
paid time scheduled, total agent preference, etc. The formulation step also includes the 
expression of contact center constraints, and agent and skill requirements in all periods to 
be scheduled as mathematical equalities or inequalities. Finally, constraints that restrict 
the values of all decision variables to nonnegative values, and some variables to 

15 nonnegative integer values only are added to the formulation. 

The MDLP model of the present invention may have one or more solutions (also 
called "feasible solutions") that satisfy all contact center constraints and requirements 
included in it. The optimal solution of an MILP model is merely finding the best feasible 

20 solution in terms of its objective Sanction value among all feasible solutions. The best 
feasible solution with the most favorable objective function value is the optimal solution 
and satisfies the optimality condition. Moreover, every agent schedule satisfying the 
optimality condition in a specific workforce environment, when all contact center 
constraints and requirements in each planning period are satisfied, is also a global or 

25 absolute optimal solution (e.g. see Taha, 1987). The optimal agent schedule is then 
reported using the optimal solution. 

To find a global optimal solution to the MILP model of an agent scheduling 
environment, the present invention applies an enhanced optimization algorithm 
30 combining the standard technique for solving MILP models optimally, the B&C 

algorithm (Wolsey, 1998), with an RA algorithm. An MILP model contains one or more 
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decision variables whose values are restricted to non-negative integer values. In the case 
of the MILP of the present invention, for example, the number of agents assigned to a 
weekly work pattern specifying their work and non-work days cann't take a fractional 
value (e.g. it can't be 4.75 agents). The B&C algorithm first solves the MILP model 
5 (called "node zero") without integrality constraints for these decision variables. This 
version of the MILP model is referred to as the Linear Programming (LP) relaxation of 
the MDLP. If the values of the decision variables in this solution also satisfy the 
integrality constraints for all decision variables, the solution is optimal. The solution 
algorithm of the present invention then terminates, and the schedule-reporting step 
10 begins. 

If the values of the decision variables in the solution of the LP relaxation at node 
zero don't satisfy all integrality constraints, then the solution algorithm of the present 
invention calls the RA algorithm to search for an integer feasible solution. The RA 

1 5 algorithm first retrieves the values of the decision variables found by the B&C algorithm 
to the LP relaxation. The RA algorithm then rounds the values of the integer restricted 
decision variables with fractional values either up or down in an attempt to find an 
integer feasible solution. After rounding the values of the decision variables, the RA 
algorithm also checks the constraints of the MILP model. If the rounded solution is 

20 violating one or more constraints of the MDLP model, values of the variables are adjusted 
to restore feasibility including addition of new agents to meet the agent requirements in 
every period, and for all skill types. If no solution satisfying all constraints can be found 
with the current solution, the RA algorithm terminates and control is returned back to the 
B&C algorithm. Otherwise, if an integer feasible solution to the MILP model is found, 

25 the RA algorithm calculates the objective value for that solution, and compares it with the 
best known solution to the MILP model, if one is available. If the new integer feasible 
solution found has a more favorable objective value, the new integer feasible solution 
replaces the best known solution to the MDLP. The RA algorithm then terminates, and 
returns the new integer solution found back to the B&C algorithm. 

30 
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If the best known (integer feasible) solution was changed in the RA algorithm, the 
B&C algorithm updates the best integer solution known for its own use, and forms new 
sub-problems (nodes) by adding new constraints (Wolsey, 1998). The entire process is 
then repeated; the B&C algorithm chooses a new sub-problem, and solves its LP 
5 relaxation. If the LP relaxation has a feasible solution satisfying all of the constraints but 
some of the integrality constraints for some decision variables, the B&C algorithm 
transfers control to the RA algorithm which searches for an integer feasible solution 
through rounding and adding more agents to the schedules. The solution algorithm of the 
present invention terminates when all sub-problems in the B&C algorithm are terminated 
10 (sub-problem termination conditions in the B&C algorithm are given in Wolsey, 1998 ). 
That is when the optimality condition is satisfied by an integer feasible solution to the 
MELP model. When this condition is satisfied, control is transferred to the schedule- 
reporting step. 

15 In the schedule reporting step of the present invention, an optimal agent schedule 

with daily start and stop times, work and non-work days, daily break times for every 
work day, agent allocation to contact types in every period to be scheduled are reported 
for every agent tour schedule. 

20 
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BRIEF DESCRPTION OF THE DRAWINGS 



For a more complete description of the present invention and its merits, reference 
should be made to the Description of the Preferred Embodiment. 

FIG. 1 shows the process architecture for optimal workforce scheduling with the 
method of the present invention, and its computer implementation. 

FIG. 2 shows a flowchart of the MILP model generation of the invention. 
FIG. 3 shows all possible work patterns in a week. 

FIG. 4 shows a flowchart of the solution method for the MILP model of the 
invention. 
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DESCRIPTION OF THE PREFERRED EMBODIMENT 

FIG. 1 shows process architecture for optimal workforce scheduling with the 
5 method of present invention, and its computer application. It consists of several sub- 
processes including forecasting (1), agent and skills requirements planning (2), 
scheduling environment information input (3) (e.g. center information, agent skills and 
work group availability, tour and shift parameters), MILP model generator (4), optimizer 
(5), and rostering (6) processes. The main function of the forecast generator is to make 

1 0 predictions of future contact volumes, and average handling times (service times) for the 
contact types offered in each planning period during the days in the schedule interval 
(e.g. a week). This sub-process uses statistical techniques to analyze patterns in the 
historical data, and develop forecasts. Agent and skill requirements planning module (2) 
combines contact volume and average handling time forecasts with service level targets 

1 5 (e.g. 80% of incoming calls answered within 20 seconds of their receipt) to determine 
agent and skill requirements in each planning period to be scheduled. In a multi-skilled 
environment, this module determines the agent and skills mix from the available agent 
skill groups that are needed in each planning period to deliver the targeted service levels 
to all contact types. This process also takes the economies of multi-skilled scheduling 

20 environment into consideration. Contact center managers input the environment 

information including the objective criterion to be optimized, and various contact center 
parameters (e.g. number, location, operating hours), agent skill and work groups, and tour 
and shin parameters to module (3). The next two modules, the MILP model generator 
(4) and optimizer (5), use the method of present invention to develop optimal agent 

25 schedules to meet agent and skills requirements developed in module (2) using the 
scheduling environment parameters provided in module (3). The present invention 
provides a method for the MILP model generation for a workforce scheduling 
environment, and an algorithm for a near-optimal or optimal solution of the MILP model 
generated by the MILP model generator to develop agent schedules. The method of the 

30 present invention covers both non-skills based, and skills based agent scheduling 
environments. 
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The scheduling environment may involve one or more contact centers. In the case 
of multiple contact centers, each contact center may be scheduled independently to meet 
their allocated agent and skills requirements with their own agent skill groups, agent 
work groups, agent availability, operating hours, and tour and shift scheduling rules. In 
this case, the method of the invention formulates the scheduling environment at each 
contact center as an MILP model and solves it optimally to obtain an optimal schedule. 
AJtematively, a plurality of contact centers may be scheduled as a virtual contact center. 
In this case, agent and skills requirements of different contact centers are combined using 
a common time zone, and the method of invention formulates one MDLP model for these 
centers to include combined agent and skills requirements with the centers' agent skill 
groups, agent work groups, agent availability, operating hours, and tour and shift 
scheduling rules. When solved together, the optimal schedule provides a combined 
schedule for agents at all of the centers included. The unique operating characteristics of 
each contact center such as operating hours, tours and shift scheduling rules, agent skill 
and work group availability, and agent availability are satisfied by modeling them in one 
common time zone, and identifying them by contact center specific identification 
numbers. 

FIG. 2 shows the basic steps of the MELP model generation method of the present 
invention. The following description, referring to FIG. 2, is presented to explain the 
MELP model generation method in its general case, before proceeding to the more 
specific scheduling environments. 

Step 7: Acquire agent and skills requirements for each planning period to be scheduled. 

Step 8: For each contact center, obtain daily operating hours, agent skill groups, 

work groups for these agent skill groups specifying various work patterns such as 
full-time 40 hours/week, or part-time 30 hours/week, and the scheduling paradigm 
as to whether to schedule each contact center separately or a plurality of them in a 
virtual contact center environment. 

Step 9: Acquire tour and shift parameters and scheduling rules, and their availability to 
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agent skill and work groups at different contact centers. 
Step 10: Acquire the objective criterion to be optimized and associated cost and benefit 
parameters. 

Step 1 1 : Generate the appropriate MGLP model by first defining decision variables, 

parameters and sets, and then formulating the objective function, agent and skill 
requirement constraints, break constraints, days-off and work pattern constraints, 
agent availability constraints, other schedule related constraints, and non- 
negativity and integrality constraints. 

Step 12: Input the variables, objective function, and constraints of the MELP model 
forniulated to optimizer (2). 

A number of MDLP formulations of the present invention for agent scheduling in a 
non-skills based (i.e. single skill) environment and their extensions are discussed first._ The 
MDLP model generation method of the present invention for skills based agent scheduling 
environments involving agents with different skills serving multiple contact types is 
disclosed after this discussion. The method of the invention for optimal solution of the MILP 
model is disclosed thereafter. 

To facilitate the description of the MILP model and its extensions used by the 
present invention in this patent document, any variable or set defined with a subscript is 
shown as flat text when the variable or set itself is a subscript of another term. For 

instance, a set defined as Bliui, is typed as BikUi in ZteBiidh where it is a subscript of the 

summation sign S. Moreover, symbols "e" and "g" are both used to imply membership 
in a set or indicate an element of a set. 
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Workforce Scheduling with Tours Requiring Consistent 
Daily Start Time and Shift Length in 
a Single Skill Environment 

5 To facilitate the presentation, this discussion assumes that all agents have the 

same skill set (i.e. not a skills based scheduling environment). In addition, the following 
assumptions regarding the scheduling environment are made to facilitate the presentation: 
Each agent works 5 days and receives two consecutive days-off during a week; There are 
one or more contact centers operating 24 hours or less a day and 7 days a week. A day is 

10 divided into equal planning periods (e.g. 15-minute intervals) that are indexed 
sequentially in Th for day h, h = 1, 7, As disclosed below, none of these 
specifications is restrictive on the MILP model of the method. A person with ordinary 
skills at art can easily extend the MDLP model of the method to cases where the number 
of work days may be less or greater than 5 days per week, a work tour less than or 

1 5 exceeding 7 days, any number of consecutive days-ofF, and a scheduling environment 
where the time period may be less or greater than 15 minutes. 

Moreover, for this discussion, it is also assumed that a tour specifies a weekly work 
pattern that is seven days long with the following breaks in the work schedule: (i) Daily 

20 breaks in the shift schedule including two relief breaks, and one lunch break during a 
work day, and (ii) two consecutive non-work days (i.e. days-off). We assume that each 
relief break is one planning period, and lunch break is two planning periods long. 
Extensions of the MELP model of the method to tours with more or less breaks,* and 
different break durations are discussed with other cases below. Agents assigned to a tour 

25 have the identical daily shift start times and shift lengths during a week. The present 
invention extends the MBLP model to the scheduling environments involving non- 
identical daily shift start times and shift lengths later in this document. The scheduling 
environment may include more than one tour type. Each tour type specifies a set of daily 
shift start times, and a shift length for seven days. In the general case, different tour types 

30 may also have different number of consecutive or non-consecutive off^ days, number of 
daily breaks and break durations. 
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Let K be the set of all distinct tour types available to assign agents. Each tour 
type k has a predetermined set of start times given in Ik. . Let aidht be equal to one if - 
period t on day h is in the shift span (that is, a work, break, or days-oflF period) of tour k 
with a daily start time of i, and zero otherwise. Suppose that the maximum number of 
agents available to assign to tour k is limited to Dr"'''''. Let Xid be the number of agents 
assigned to tour k with a daily start time of period i. Let Cki be the weekly cost (wages 
and benefits, etc.) of assigning an agent to tour k with a daily start time of i, ielk , and 
ksK. Let Ukiin, Wkiht and Vkiht, daily break variables, be the numbers of agents assigned 
to tour k with a daily start time of i and starting their first relief, lunch and second relief 
breaks in period t on day h, respectively. Daily relief and lunch breaks must start and end 
within the break windows specified for them. Let B Ikiii, B2kih and B3idii be the sets of 
planning periods on day h during which an employee assigned to tour k with a daily start 
time of i may start, respectively, his/her first relief, lunch and second relief break and 
complete within the time windows specified for these breaks. Break variables Ukih(ti), 
Widh(t2) and Vkih(G) are defined for tlEBlkUi, t2eB2kiii, t3EB3kii„ and islk, keK, h = 1,...,7 
respectively. Further, let Tlkhi, T2kht and T3ki,t be the sets of start times for tour k for 
which period t on day h is, respectively, a first relief, lunch, and second relief start period 
within the break windows specified for them. Let Ykih be the number of agents assigned 
to tour k with a daily start time of i who are scheduled start their consecutive non-work 
days (days-off) on day h. Finally, let bht be the number of agents required in period t on 
day h to deliver the desired service level to customers. The JVULP model of the present 
invention for the scheduling environment assumed can be stated as follows. 
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Xki (1) 

Subject to 

SiceK ^^lElk akihl Xfej - SkeK ^ielk akiht (Ykih + Yici(h-1)) 

- SkeK SiETIkltf Ukilu " ^keK ^isT2Vhl (Wkiht + Wkih(|.l)) 

ZkeK 2ieT3kht Vkiht > bht , 

teTh, h= 1,...,7, 
(2) 
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Xki - Ykih - Yid(h-i) = SteBikih Ukiht , ielk, keK, h = 1,...,7, (3) 



Xki - Ykih - Yki(h-i) = 2teB2kih Wkflit ,islk, keK,h=l,...,7, (4) 



Xki - Ykih - Yki(h.i) = 2teB3kiii Vkiht , islk, ksK, h- 1,...,7, (5) 



20 Xki = 2h Ykih , ielk, keK, (6) 



SieikXki < Dk""" ,keK, (7) 



Xki, Ykih, Ukiht, Wkiiu, and Vkiht ^0, and integer for all k, i, h, and t. (8) 
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Objective function (1) is the total cost of assigning agents to various tour types and start 
times per week. Constraint (2) ensures that the total number of agents assigned to various 
tours and available in period t on day h (first term on the left side of (2)) minus the 
number of agents who are assigned to these tours but taking this day off (second term on 
5 the left side of (2)) and the number of agents who are at work on day h but taking their 
first relief (third term on the left side of (2)), lunch (forth term on the left side of (2)), 
and second relief (fifth term on the left side of (2)) in period t on day h is sufficient to 
meet the agent requirement, bht, in that period. The second term on the left side of (2) 
includes agents who are starting their consecutive days-off" schedule on day h as well as 

10 those who started a day eariier on day (h-1). Break constraints (3-5) ensure that one first 
relief, one lunch, and one second relief break are scheduled for each agent who is 
assigned to tour k with a daily start time of i and working on day h. Thus, the number of 
first relief, lunch, and second relief breaks scheduled for tour k with a daily start time of i 
on a given day is equal to the number of agents assigned to this tour and start time minus 

1 5 the number of them scheduled to take day h off, Xki - Ykih - Yiu(h.i) (because of two 
consecutive days-ofif per week requirement). Constraint (6) ensures that a set of two 
consecutive days-off is scheduled for every agent assigned to tour k with a daily start 
time of i. Constraint (7) ensures that the total number of agents assigned to tour type k 
doesn't exceed the maximum number of agents that are available to assign to this tour 

20 type. This constraint can easily be modified to limit the number of agents that can be 
assigned^ a set of tour types (e.g. fijll time 40-hour a week tours may include 8-hour 
shifts and 5 work days a week, and 10-hour shifts and 4 work days a week). Finally, 
constraint (8) restricts the tour variables Xid, days-off* variables Yidh, and break variables 
Uidiii, Widin, and Vkiiu to nonnegative integers to avoid negative values and/or fractional 

25 agent schedules. 

Extensions of MILP Model (1-8) 

30 It should be clear to those skilled in the art that the specific embodiments and the 

MILP model disclosed above may be readily utilized as a basis and extended to other 
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scheduling environnients without departing from the scope of the invention. For 
example, objective function (1) can be easily reformulated to minimize the total amount 
of agent time scheduled or paid agent time scheduled by setting Cid equal to scheduled or 
paid time per week per agent assigned to tour k with a daily start time of i, islk, and keK. 
5 Also, if agents rate their most preferred, second most preferred, etc.', tours using an 

RAting system, these ratings for various tours and start times can be aggregated and used 
in the objective function to maximize agent preference in scheduling tours and start 
times. Note that in this case, in order to not to schedule redundant tours (a redundant 
tour is a tour that can be removed without causing agent shortages in any of the periods 
10 covered by that tour during the week) and increase agent costs unnecessarily, rating 

system should assign the lowest value to the most preferred tour, a higher to the second 
most preferred tour, etc. Consequently, the objective function will remain as of 
minimization type to maximize agent preferences. 

15 Another extension of MCLP model (1 -8) involves agent scheduling when there 

aren't enough agents to meet the agent requirements in every period to be scheduled. In 
this case, objective fiinction (1), and constraint (2) are modified by adding overstaffing 
and understaflfing variables. Let Ohtbe the number of agents in excess of agent 
requirements, and Sht be the agent shortage in period t on day h. A penalty for shortages 

20 in a period, Pht, is included in the objective function to find the optimal schedule that 
minimizes the total agent cost plus the shortage cost. 

Minimize Zk^K ^itik Cki Xki + Zh Eieji, Pht Sht (9) 

25 Constraint (2) is then formulated by adding Om for overstaffing and Sm for shortages to 
its lefl side as follows; 

ZkeK Zieik akiht Xki - ZkeK ZiEik akiht (Ykih YkiXh-l)) 
30 - SkeK ZieTlkht Ukiht " ^keK ZieT2kht (Wkiht + Wkih(t-1)) 
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- 2)keK 2]iET3kht Vkiht + Shi - Oht = bht , 

tsTh, h= 1,...,7, 

(10) 

Sht , Oht ^0, for all tsTh and h. 

If a certain amount of overstaffing is needed for other tasks (e.g. to answer email 
or fax based contacts), it can be ensured by restricting the overstaffing to be greater than 
or equal to the desired level fht. 

^tsTh Oht > fht , tsTh, h= 1,. . .,7, (1 la) 

If up to a certain amount of understafEng staffing is tolerable (e.g. to see if the 
total cost can be reduced), this can be achieved by adding a similar constraint: 

StsTii Shi > dht , tsTh, h= 1,. . .,7, (1 lb) 

where dht is the maximum understafiBng tolerated on day h. Constraints (11a) and (lib) 
can easily be further modified to impose limits for a specific period(s) on day h. 

If a tour type has more than three breaks on a given day (e.g. first relief, lunch, 
second, and third relief breaks on Thursdays and Fridays), the method of the present 
invention extends MILP model (1-8) by introducing a set of break variables and a break 
constraint (similar to constraints (5-7)) for every break that needs to be scheduled so that 
every agent assigned to tour k and working on that day will have the correct number of 
breaks scheduled for him/her. If agents assigned to a tour are given fewer breaks (e.g. 
part time agents working 4 hours a day), the method of the present invention extends 
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MILP model (1-8) by excluding the break variables and constraints associated with the 
breaks not allowed. 

Furthermore, if a tour type has a break longer than 2 consecutive planning 
5 periods, then the method of the present invention extends MILP model (1-8) by 

subtracting all agents who started their breaks in an earlier period but haven't completed 
yet and the agents who are starting their break in period t from the left side of constraint 
(2). For a 4-period break, for example, agents who started their breaks three periods 
earlier, two periods eariier, one period earlier, and agents starting in the current period 
1 0 will be subtracted from the number of agents scheduled in that period from the left side of 
constraint (2). 

In MILP model (1-8), even though daily shift start times for an agent remain the 
same on different days, break times may vary from one work day to another. For 
1 5 instance, it is possible to have the first relief break of an agent be scheduled at 9:00am on 
Monday, and at 10:00am on Tuesday. If a tour (say tour k) requires consistent relief and 
lunch break times on all work days in addition to consistent daily shift start times and 
shift lengths, these requirements are ensured by adding the following constraints the 
MILP model (1-8): 

20 

2h Uidjit £ qk SUkit , 
: ShWidht < qkSWkit, 

2h Vkiht < qk SVkit , 

25 2iteDlki SUkii 1, 

2tED2ki SWkit = 1, 
StED3kiSVkii =1, 

SUkit, SWkit, SVkit > 0 and integer. 



,isIk,tsDlki, (I2a) 
, islk, tsD2ki , (12b) 
, ielk, tsD3ki , (12c) 



,isIk,kGK, (13a) 
,ieIk,kGK, (13b) 
, islk, kGK, (13c) 
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where Dlid , D2ki , and D3ki are the sets containing the first consistent (across the days 
of a week) daily first relief, lunch, and second break times in Blkih, B2kiii, and B3kih, 
respectively, qk is a very large number, SUkii, SWidt, and SVi^t are binary variables taking 
5 a value of either zero or one to indicate whether period t is selected as the consistent daily 
first relief, lunch, or second relief break time, respectively. 

Constraints (12a) and (13a) are restricting the daily first relief break times to one 
of the break times specified in Blidii. Constraint (13a) ensures that only one relief break 

10 time is selected, that is, only one SUidt is equal to one while the others are zero. Thus, for 
example, if the first relief break times for a tour starting at 7:00am daily are BlkUi = 
{8:15am, 8:30am, 8:45am}, h= 1,...,7, then only one of SUki(8:i5) , SUici(8:30) , or SUki(8:45) 
will be equal to one and the others will be equal to zero. Likewise, constraints (12b) and 
(13b) are restricting the lunch break times, and constraints (12c) and (13c) are restricting 

_15 .the second. break times-to. consistent dailyschedules. . . - — 

MILP model (1-8) of the method is presented for tour types that require two 
consecutive non-work days be scheduled for each agent. The present invention extends 
MILP model (1-8) to include tours with 3 or more consecutive days-off by subtracting 
20 fi-om the left sides of constraints (2-5) the agents assigned to this type of tours and who 
are off as part of their scheduled consecutive non-work days, and agents starting their 
days-oflf on day h. 

If, however, tours in K require a minimum of two consecutive days off and, in 
25 addition, a total of rk non-consecutive days off, the present invention extends MLLP 
model (1-8) as follows: 



30 
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Minimize SkeK ^ieik Cki Xkj (14) 
Subjecto 

^keK 2^ielk aidht Xki - ^keK ^'isJk Skiht (Yjah + Yki(h-i)) 

- SkeK ^islk 2mf=lu(h.l) Skiht Zkimh " ^keK ^isTlkht Ukihl 

- SkeK 2]i£T2klit (Wkiht + Wkih(M)) " ^^keK 2ieT3kht Vfciht > bht 

,t8Th,h=l,...,7, (15) 
Xki - Ykih -Yki(h-l) - 2^m=^h,(Ii-l) Zkimh = SteBlkih Ukiht ,isIk,keK,h=l,...,7, (16) 



Xici - Ykih -Ykj(h-i) - 2m*h,(ii-i) Zkimh = ^eB2kih Wkiht. , is^c, keK, h = I,. -, 7, (17) 

Xki - Ykih -Yki(h-1) - 2m=*h,(h-l) Zkimh = 2teB3kili Vkiht , islk, keK, h = 1,...,7, (18) 

Tk Ykih = 2i^(h+i) Zkihi , islk, keK, h = 1,...,7, (19) 

Zkihm < Ykih , islk, keK, h = 1,...,7, 

mf=h, (h+1) (20) 

Xki = Eh Ywh , ielk, keK, (21) 

Si^ncXki < Dk"^ , keK, (22) 

Xki, Ykih, Zkiimi, Ukiht, Wkiht, and Vkiht >.0 and integer for all k, i, h, m, and t. (23) 
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where Zkihm is defined as the number of agents who are assigned to tour k with a daily 
start time of i and will start their consecutive days off on day h and will also take day m 
off as part of their tour schedule. Note that Zidimi is defined for m =# h, (h+1), for tours 
requiring a minimum of two consecutive days off, in order to avoid infeasible days off 
5 schedules. 

Constraint (19) ensures that each agent assigned to tour k receives exactly rk 
additional days off. If for tour k, a maximum of (2 + r^) days-oflf with a minimum of 2 
consecutive days-oflf, this can be formulated by simply modifying constraint (19) as a 
10 greater than or equal to type constraint (that is, the term on the left side of (19) is greater 
than or equal to the term on the right side). 

Constraints (19) and (20) ensure that r^ different days are scheduled as non-work 
days for every agent on days other than their consecutive days off. In constraints (15-18), 
1 5 agents who are oflF on day h due to their consecutive days-off schedule as well as the 
agents who are oflF due to their additional days off schedule are subtracted from the total 
number of agents scheduled. 

When a scheduling environment involves only tour types not requiring any 
20 consecutive days off, that is equivalent to requiring a minimum of one "consecutive" day 
off. In this case, the present invention extends MELP model by removing Zjdhi variables 
from all constraints, and dropping constraints (19) and (20). Ykih is now defined as the 
number of agents assigned to tour k with a daily start time of i who will take day h off. 
Note that only Yidi, will be in constraints (15-18) instead of both Ykii, and Yid(h-i) since the 
25 tour types considered do not require consecutive days off. The MELP model is now 
stated as follows: 



30 



25 



Best Available Copy 

Minimize ZkeK ^\ejk Cki (24) 
Subjec to 

5 SkeK SiEDc akiht ^ " 2keK ^islk akiht Ykfl, - SfceK ^igTlkht Ukiht 

- SkeK 2isT2kht (Wkiht + Wkih(t-1)) " ^keK 2isT3kht Vkiht > bhi 

,t8Th,h-l,...,7, (25) 

10 

Xki - Ykih = SteBikih Ukiht , islk, keK, h = 1,.,.,7, (26) 

Xki - Ykih = 2,eB2kih Wkiht , islk, keK, h = 1,...,7, (27) 



15 Xki - Ykih = S,eB3kih Vkiht , ielk, keK, h = 1,...,7, (28) 

Fk Xki = 2h Ykih , ielk, keK, (29) 

2i,ik Xki <. Dk"""" ,keK, (30) 

20 

Xki, Ykih, Ukiht, Wkiht, and Vkiht >.0 and integer for all k, i, h, and t. (31) 
where rk is now defined as the number of days-oflF an agent assigned to tour k receives. 

25 

A tour is said to have a fixed day off if every agent assigned that tour will be off 
on the same day. For example, a contact center may close on Sundays. In this case, all 
agents will receive a day off on Sunday. The present invention extends the MILP model 
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by setting the days-ofF variables for the fixed days off equal to the number of agent 
assigned to that tour. For example, if tour k requires two consecutive days off, and day h 
is an off day for all agents assigned to this tour, then 

Xid = Yidh + Yki(h-i) jElk, (32) 

And if tour k doesn't require consecutive days off, and day h is an off day for all agents 
assigned to this tour, then 

Xki = Ykih ,islk. (33) 

If operating hours are shorter for certain days (e.g. regular hours of 5:00am to 
10:00pm with Saturday operating hours of 5:00am to 4:00pm) for a tour or a contact 
center, the present invention extends the MILP model by formulating each such tour as 
consisting of two tours with identi cal w ork rule s exce pt for their start times as follows: (i) 
tour 1 includes the start times for which the latest daily shift completion time not 
exceeding the operating hours for the early closure day, and (ii) tour 2 has the start times 
not covered by tour group 1 (ending after closure on the early closure day) and has a 
fixed day off on the early closure day to not to schedule any agent to work on those days, 

A scheduling environment can have any combination of the extensions and tour 
types discussed above. The MDLP model of the present invention can be extended to 
formulate these scheduling environments by combining various extensions of the 
constraints and objective fijnction together with the definitions of decision variables and 
parameters disclosed above. 
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Workforce Scheduling with Tours not Requiring Consistent 
Daily Start Time and Shift Length in 
a Single Skill Environment 



To facilitate the discussion and notation, the scheduling environment described 
before the disclosure of MDLP model (1-8) is considered again. The present invention 
also extends the MELP model to tours not requiring consistent daily shift start times and 
same daily shift length requirements. For these tours, shift start time as well as shift 
10 duration for an agent may vary fi-om one day to another. 

Assume that the scheduling environment involves only tours not requiring 
consistent daily start time and shift length. Let QK be the set of all such tour types. 
Assume that agents assigned to tour type k can start at any one of a set of predetermined 
1 5 start times given in Qlk. Let the daily shift lengths allowed for tour type k be Fk. Fk may 



also change from one day to another. Note that tour k has a minimum weekly work limit 
specified by the. shortest shift in Fk, and a maximum weekly work limit specified by the 
longest shift in Fk. The daily shift start times should not allow for back-to-back shift 
schedules for agents without allowing enough time for rest. Thus, the latest start time for 
20 a tour plus the length of the longest daily shift type are not allowed to exceed 24 hours 
minus the desired minimum rest period between consecutive daily shifts. To prevent this 
for a tour type, the present invention creates separate tour types each v^th start times 
allowing a minimum required time between the latest possible end of a shift on one day 
and the earliest possible start of a shift on the following day. 

25 

Moreover, assume that the maximum number of agents available to assign to tour 
k is QDk"™^. Let QX^ be the number of agents assigned to tour k and shift length n 
with a starting of period i on day h. Define daily break variables QUbiiht, QWkniiit and 
QVknihi as the numbers of agents assigned to tour k and shift length n with a daily start 
30 time of i and starting their first relief, lunch and second relief breaks in period t on day h, 
respectively. Depending on the shift length, some of the breaks may not be provided. 
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This case is formulated in the MILP model as disclosed before for tours with more/fewer 
daily breaks. Daily relief and lunch breaks must start and end within the break windows 
specified for them. Let QBlkmii, QB2kmh and QBSMibe the sets of planning periods on 
day h during which an employee assigned to tour k and shift length n with a daily start 
5 time of i may start his/her first relief, lunch and second relief break and complete within 
the time windows specified for these breaks, respectively. Break variables QUknih(ti), 
QWk„iiKt2) and QViouiKO) are defined for tieQBlknih, t2GQB2kmh, and t3GQB3iaiih, 
respectively. Further, let QTlknht, QT2kri,t and QTSiojit be the sets of start times for shift 
length n for tour k for which period t on day h is a first relief, lunch, and second relief 

10 start period, respectively. Let aicniht be equal to one if period t on day h is in the shift span 
(that is, a work or a break period) of agents assigned to tour k who have the start time of 
i and shift length n on day h, and zero otherwise. A work pattern in this case is a string 
of zeros and ones indicating work days (=1) and non-work days (=0) scheduled for agents 
assigned to this work pattern. All possible work patterns for a week are shown in FIG. 3 

15 Let the set of all allowed w ork and non-work patterns for the ag ents a s signed to tour k be 
QLk. Only the work patterns in FIG. 3 satisfying the work and non-work day rules 
specified for tour k are included in QLk (e.g. a minimum of two consecutive days off, no 
work on Sunday, etc.). A work pattern in this case specifies work and non-work days 
explicitly. Define the days-off pattern variable Qki as the number of agents assigned to 

20 tour k with a work pattern of 1. Let Akjh be equal to one if day h is a work day for 
agents assigned to work pattern 1 of tour k, and zero otherwise. 

Let Ckj be the weekly "fixed" cost (e.g. benefits, etc.) of assigning an agent to tour k v^th 
work pattern of 1, leQLk , and ksQK. Let Ckni be the daily wage paid to agents assigned 
to tour type k with a daily shift length of n and start time of i. Finally, let bht be the 
25 number of agents required in period t on day h to deliver the desired service level to 
customers. The present invention extends the MILP model to this scheduling 
environment as follows. 
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Minimize ZkepK ^)sQlk Cu Qkl + SkeQK SnePk SisQDc Ckni QXkrfj (34) 

Subjec to 

SkeQK SnsFk SisQlk aknilit QXknhi " SkeQK SnsPk ^iEQTIknht QUkniht 

- SkeQK 2nsFk 2iEQT2knht (QWkniht + QWkjuH(t-l)) 



- ^^keQK SnsFk -SjEQTSknht QVknihf > bht 

10 ,tsTh,h=l,...,7, (35) 

QXknhi = SteQBiknih QUkniht , nsFk, IsQIk, keQK, h = 1,...,7, (36) 



QXkriii = EteQB2kmh QWfcmhi , nsFfc, isQIfc, ksQK, h = 1,...,7, (37) 

15 . 

QXknhi = IlteQBSknih QVfcnibt , HsFk, IsQIk, ksQK, h = 1,...,7, (38) 



SiEQuAkihQki = ^^iKFk SiEik QXkrili , keQK, h= 1,...,7, (39) 



20 IieQuQki < QDk""" , keQK, (40) 



QXkniii, Qkl, QUkniht, QWkniht, ahd QVknihi >_0 and integer for all k, n, i, h, and t. (41) 

25 

Expression (34) above is the objective fijnction that minimizes the total "fixed" 
cost of assigning agents to tours with various feasible work patterns plus the total wages 
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paid to agents. Constraint (35) ensures that the total number of agents assigned to various 
tours and shift lengths, and available in period t on day h minus the number of agents 
who are at work on day h but taking their first relief (second term on the left side of 
(35)), lunch (third term on the left side of (35)), and second relief (four term on the left 
5 side of (35)) in period t on day h is sufficient to meet the agent requirements in that 
period. Break constraints (36-38) ensure that all required breaks are scheduled for each 
agent who is assigned to tour k. These constraints as well as the break variables are 
defined only for valid breaks for a shift length. Constraint (39) ensures that the number 
of shifts scheduled on day h is equal to the number of agents who work on that day based 

10 on the scheduled work patterns. Constraint (40) restricts the total number of agents 

assigned to different feasible work patterns of tour type k to be less than or equal to the 
total number of agents available to assign to this tour. This constraint can easily be 
modified to limit the number of agents that can be assigned to a set of tours. Finally, 
constraint (41) restricts the work pattern, daily shift, and break variables to nonnegative 

1 5 integers to avoid negative values and/or fractional agent schedules. 



Extensions of MILP Model (34-41) 

20 Objective fimction (34) can be easily reformulated to minimize the total amount 

of agent time scheduled (or paid agent time) by setting Cki equal to scheduled time (or 
paid agent time) per agent assigned to tour k and work pattern 1, leQLk, and keQK. 
Similarly, if agents rate their most preferred, second most preferred, etc., tours using a 
rating system, these ratings for various tours and work patterns can be aggregated and 

25 used in the objective value to maximize agent preference in scheduling tours and work 
patterns. In order to not to schedule redundant tours and increase agent costs 
unnecessarily, rating system should assign the lowest value to the most preferred tour, 
and next value higher to the second most preferred tour, etc. Consequently, the objective 
fiinction will remain as of minimization type to the maximize agent preference. 

30 
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The MILP model can be easily modified to develop optimal schedules when there 
are not enough agents to meet the agent requirements in all periods (i.e when there are 
going to be agent shortages in some periods). In this case, objective function (34) is 
modified as follows: 

5 

Minimize ZieQlIc ^^IcsQK Ckl Qkl + ^keQK ^neFk ^ieQIk Ckni QXknhi 

+ ZteQTl. Ish Pht Shi (42) 

10 

Constraint (35) is then formulated by adding Oht for overstaffing and Sht for shortage to 
its left hand side as follows: 



15 

SkeQK S^nsFk ^isQik akiiiJit QXknhi - ^keQK ^^nsFk SieQTlknlit QUkiiihl 



2k€QK 2n£Fk SieQT2kiilit (QWkniht + QWkiiih(t-l)) 



20 - S^keQK ^neFk ^isQTSknht QVkniht + Shi - Oht ^ bht 



,t8Th,h=l,...,7, (43) 



25 Sht , Oht ^0, for all tsQTh and h. (44) 

If a certain amount of overstaffing is needed for other tasks (e.g. to answer email 
or fax based contacts), it can be assured by restricting the overstaffing to be greater than 
or equal to the desired level fju. 
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teTh, h= 1,...,7, 



(45) 



If up to a certain amount of understaffing staffing is tolerable (e.g. to see if the 
5 total cost can be reduced), this can be achieved by adding a similar constraint: 




tsTh, h= 1,...,7, 



(46) 



where dht is the minimum overtime desired on day h. Note that (45) and (46) can easily 
10 modified to impose limits for a specific period(s) on day h. 

This case is similar to one of the cases discussed earlier for formulation (1-8). 
Constraints (2-5) and break variables in (1-8) have the same structure with constraints 
(35-38) and break variables. In both cases, the number of daily breaks scheduled is equal 
1 5 to the number of agents assigned to a shift and working that day. Thus, the existence of 
non-consistent daily start times and shift lengths doesn't affect the way breaks are 
scheduled in the model. Consequently, tours with more breaks or longer breaks can be 
modeled as disclosed before for formulation (1-8). 

20 Since the valid work patterns are determined in advance in formulating (34-41), 

the minimum, maximum, and exact number of consecutive days off* and/or total number 
of days off requirements for a tour type can be easily included in fonnulation (34-41). 
All possible work patterns are shown in FIG. 3. All feasible work patterns satisfying the 
days-off rules for tour k, keQK, are then identified from this list and included in QL^. 

25 Their coefficients are included in constraint (39). 

If a tour type has fixed days oflF for all agents assigned to it then all feasible work 
patterns selected for that tour should have a non-work day on that day. To reduce the size 
of the formulation, shift and break variables for the fixed days off as well as associated 
30 break constraints (36-38) are excluded from formulation (34-41). 
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If working hours are shorter on certain days for a tour or a contact center, the 
present invention extends the MILP model by excluding the daily shift variables that has 
a late start time resulting an end time after the early closure time for that tour or center. 
Break variables as well as break constraints (34-41) associated with these shift variables 
are excluded. 



Workforce Scheduling with Tours Requiring Both Consistent 
Daily Start Times and Shift Lengths, and Non-Consistent 
Start Times and Shift lengths in 
a Single Skill Environment 



A workforce scheduling environment may include tours with and without 
consistent daily start time and shift length requirements discussed earli er. The present 
invention extends the MELP model to this type of workforce scheduling environment in 
two ways. First, even though formulation (34-41) is disclosed earlier for scheduling 
environments involving only tours not requiring consistent daily start times and shift 
length, it is extended to model such tours with a slight modification as follows; If a tour 
requires consistent daily start times and shift lengths, then each daily start time for this 
tour type represents a pseudo-tour type that has the same start time and shift length on all 
days. Thus, a tour that can start every 30 minutes in a 7-day-a-week-and-24-hour-a-day 
operation will have 48 (for 48 daily start time) pseudo tour types. Each pseudo tour k has 
only one daily start time in the daily start time set QI^. Consequently, formulation (34- 
41) will have work pattern variables Qu, and one daily shift variable QXkhi, and the 
associated daily break variables for that shifts for tour k, and day h = 1,...,7, for each 
such pseudo tour. 

Moreover, if a tour has multiple shift lengths but doesn't allow varying daily shift 
start times, the present invention extends the MELP model by introducing pseudo tour 
types (one per daily start type) and formulating the shift and work pattern variables as in 
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formulation (34-41) with daily shift lengths in Fk. This way both consistent daily start 
time requirement and aviailability of multiple daily shift lengths are included in the MILP 
model. 

Another way that the present invention extends the MILP model is to merge 
formulations (1-8) and (34-41) for this type of scheduling environments to obtain optimal 
schedules. To distinguish cost parameters, Cx:*ki is introduced as the cost of assigning an 
agent to Xid . Using the definitions of the parameters and variables given for 
formulations (1-8) and (34-41), merging the objective functions (1) and (34) and demand 
constraints (2) and (35), and including the necessary break and days-off constraints from 
formulations (1-8) and (34-41), the present invention extends the MELP model for a 
scheduling environment involving both consistent and non-consistent shift start times and 
shift lengths as follows: ' 

-Minimize— Zk6K-2isik--Cxid-Xki =^—2]keQK-2ieQUc-Cid-Qki — 



'keQK ^ncFk ^isQDc Cjcni QXicnhi 



(47) 



Subject to 




- ZkeK 2isT2kht (Wkiht + Wkih(M)) " ^keK ^ieTSkht Vidht 



+ 2, 



fkeQK ^iicFk ^icQIk akiiihl QX|aiiii 



- Xi 



'keQK ^nsFk^ieQTlknhl QUkiiiht 



-2, 



'keQK SnsFk 2iEQT2knht (QWkniht + QWknih(t-l)) 




t8Th,h= 1,...,7 



(48) 
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Xki - Ykih - Yki(h.i) - SteBikih Ukiht , islk, ksK, h - \,...J, 



Xki - Ykih - Yki(h-i) - 2,6B2kih Wkihi , ielk, keK, h - 1,...,7, 



5 Xki - Ykih - Yki(h-i) - 2teB3kih Vfciht , isIk, keK, h - 1,...,7, 



Xki = Ykih , ielk, keK, 



10 



QXknhi = litgQBikmh QUkniht , nsFk, ieQIk, keQK, h = 1,...,7, 



QXkniii = I)teQB2kmh QWkiuitt , nsFk, isQIk, keQK, h= 1,...,7, 



QXkriii = 2tsQB3kiiih QVkniht , neFk, isQIk, keQK, h= 1,...,7, 

15 SiEQLk Akih Qki = SnsFk SisDc QXfcnhi , keQK, h = 1,...,7, 

SicikXki < Dk"'" ,keK, 

IisQuQki < QDk™'^ , keQK, 

20 

Xki, Ykih, QXkriii, Qki, Ukiitt, Wkiht, and Vkiht, QUkniht, QWkmht, and QVkniht >_0 and 
integer for all k, n, i, h, and t. 
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where objective function (47) combines objectives functions (1) and (34), constraint (48) 
combines agent, break, and days-ofF variables specifying agent availability in a period 
from constraints (2) and (35), constraints (49-51) are the break balance constraints for 
the tours requiring consistent daily shift start times and shift lengths, constraint (52) 
5 ensures that sufficient number of days-off are scheduled for the agents assigned to tours 
requiring consistent daily shift start times and shift lengths, constraints (53-55) ensure 
that sufficient number of first, lunch, and second relief breaks are scheduled for the 
agents assigned to tours not requiring consistent daily shift start times and shift lengths, 
constraint (56) ensures that the number of shifts scheduled on a given day is equal to the 

1 0 number agents who are assigned to work on that day by the work patterns scheduled, 
constraints (57) and (58) ensure that the numbers of agents assigned to various tours 
requiring consistent or non-consistent daily shift start'times and shift lengths don't exceed 
the limits on the agent availability for these tours. Finally, constraint (59) restricts the 
tour, days-ofF, daily shift, work pattern, and break variables to nonnegative integers to 

1 5 avoid negative values and/or fractional agent schedules. 



Depending on the scheduling environment and its tour and shift requirements, the 
present invention extends the MILP model by merging formulations (1-8), (14-23), (24- 
31), (34-41) or any extensions of these as disclosed before in a manner similar to 
formulation (47-59) above. It should be clear to those skilled in the art that the specific 
embodiments and the MELP model disclosed above can be readily utilized as a basis and 
extended to other scheduling environments involving the types of tours considered in 
formulations (1-8), (14-23), (24-31), and (34-41), and their variations without departing 
from the scope of the invention. 

Workforce Scheduling in a Skills Based Environment 

The present invention extends the MELP model to scheduling environments 
involving multiple agent skill groups, agent work groups, and multiple contact types 
requiring different agent skills. An agent skill group contains agents with the same skill 
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sets. Let R be the set of different contact types with known agent requirements in each 
period t and day h. Let the number of agents with qualified skills for contact type reR 
required in period t on day h be b^hi- Suppose now that there are a number of agent skill 
groups jsJ, each with a unique set of skills. Agent requirements should reflect the multi- 
skill efficiency expected in a skills-based routing environment due to agents with 
multiple skills and, depending on the scheduling environment, may be specified for skills 
groups rather than contact types. Let the set of contact types that skills group j is 
qualified to provide service be Nj, jeJ. A contact type may be served by, multiple agent 
groups. Let the set of skills groups that can service contact type r be M". 

The present invention extends the MILP model for skills based scheduling by 
formulating the agent assignment to contact types, and contact type based agent 
requirements fulfillment in separate constraints. To facilitate the description, consider 
formulation (34-41) and constraint (35). Further, consider the availability of agents from 

group j . The number of a g ents that are available from group, jjn^p_eriodjjo>nday_h_is 

equal to the number of agents scheduled to work on day h minus the number of agents 
taking a relief or lunch break in period t. Let G'^'ht be the number of agents fi-om skill 
group j allocated to serve type r contacts, where reNj, and jeM^ These variables are also 
referred to as "skill allocation variables". Thus, 

2^k€QKj 2neFk 2]ieQIk akniht QX'iaihi " SkeQKj ^neFk^^igQTlknht QU'kniht 

(Qww + QwWd) 

- ^keQKj ^neFk ^iEQT3knht QV'kniht " ^reNj - 0> 

j6J,t8Th,h=l,...,7, (61) 

The first term on the left side of (61) is the number of agents fi"om group j scheduled in 
period t on day h. The second through fourth terms on the left side are the numbers of 
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agents who are scheduled to be on a relief or lunch break during this period. Thus, the 
difference is equal to the number of agents who are scheduled and working (available) 
during period t on day h. Finally, the fifth term is the sum the numbers of agents 
allocated to various contact types in Nj, 

The present invention extends the MILP model into the skills based scheduling 
environment in two different ways. First approach involves merging formulations (14- 
23) and (34-41). Consider constraints (15) and (35). In the skills based scheduling 
environment, the right sides of these two constraints are merged as follows: 

f jht(X, Y, Z, U, W, V, QX, Q, QU, QW, QV) = Z^eKj Sieik a^iht X^id 

- StceKj SiElk aidht (Y-*kih + Y^ki(h-1)) - ^keKj ^^isDc 2^m^(h-l) ajdht Z'ldmh 
=-2keKj-2isT4khl-U!kmt^_2keKj-2ieT2kht-(^l^^ 

- ^keKj 2isT3kIit Vkiht 

+ SkeQKj 2]neFk ^isQDc akniht QX-'knhi 

- ^keQKj 2^ii£Fk ^^isQTlkiilil QU^kniht 

- I^keQKj SneFk ^ieQTlknht (QW^knihl QWknih(M)) 

- ^^keQKj SneFk 2igQT3knht QVkniht 

j8J,tsTh,h=l,...,7, (62) 

The terms on the right hand side of (62) represented by fhi(X, Y, Z, U, W, V, QX, Q, QU, 
QW, QV) give the number of agents from group j assigned to various tours (with and 
without consistent start time and length) and are available in period t on day h. The 
notation introduced for (1-61) is now extended to distinguish agent skills groups. For 
example, X = { X^ki: ielk, , keKj, h = 1,...,7} where X-'ki is as defined before except it now 
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represents the number of agents from group j assigned to tour k with a daily start time of 
i. The sets Y, Z, U, W, V, QX, Q, QU, QW, and QV are defined similar to the set X to 
include the respective variables for the agent groups in J, Constraint (61) is now 
extended to the skills based scheduling MILP model of the present invention as follows. 

fVX, Y,Z,U,W, V,QX,Q,QU,QW,QV)- Z,eNj (^\t = 0, 

,j= 1, ...,n,t8Th, h- 1,...,7, (63) 

Given the agent allocation G^'^hi in every period and day for different contact 
types, agent requirement constraints for contact type r are formulated in the MILP model 
as follows. 

SjsMr ei'G^m > b^hi, ,rER,teTh,h= l,...,7, (64) 



where e*^ is the relative efficiency of an agent from agent group j in serving contact type 
r with respect to an agent whose primary skill (highest proficiency level) is serving 
contact type r, e'^8[0, 1]. For instance, if e**^ = 1, agent group j is equally efficient, and if 
d*^ = 0.75, agents in group j are 75% efficient with respect to an agent whose primary skill 
is serving contact type r. Values of G'^hi variables do not have to be integer. 

Variable definitions presented before are modified to include agent group jeJ as a 
superscript in all variables. To schedule breaks and days off for different agent groups, 
constraints (16-18) and (36-38), except for variable definitions for different agent skill 
groups, remain the same since the agent allocation does not affect the break and days-off 
scheduling constraints. Objective functions (14) and (34) are merged and modified to 
include potential wage differentials for different agent groups. The MELP model of the 
present invention is now disclosed below. 
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Minimize ZjeJ Sislk SkeKj Cj^U X'ki + ^^jeJ^^lsQLk SkeQKj Ckl Q^kl 

+ SjtjSkeQKj SiisFk ^ieQIk c'kni QX'knhi (65) 

5 Subject to 

SjsMr > b'h. ,reR,t8Th,h=l,...,7, (66) 

f Jht(X, Y, Z, U, W, V, QX, Q, QU, QW, QV) 

JsJ,teTh,h=l,...,7, (67) 
X'ki - Y'kih - Y'ki(h-i) - I^m=>h.(h-i) Z W = ZisBikih Lfkiht , jeJ, ielk, keKj h = K -J, (68) 



15 X'ki - Y'kih - Y^ki(h.i) - 2m#h,(h-i) Z^kimh = 2teB2kih Wfciiit , jsJ, ish, keKj, h = 1,...,7, (69) 
X'ki - Y-^kih - Y-'kj(h-i) - 2m*h,(h-i) Z^kimh = 2t6B3kih Vfcjht , jsJ, Islfc, ksKj, h = 1,...,7, (70) 
Tk Y^kih = 2iHi(h+i) Z^kihi JeJ, ielk, keKj, h = 1,...,7, (71) 

20 

< Y^kih , jeJ, ielk, ksKj, h = 1,...,7, 

m=/=h,(h+l) (72) 
X^ki = 2h Y\ih JeJ, ielk, keKj, (73) 

25 QX^knhi = 2. t€QBlknih QLPknihi , jeJ, neFk, ieQIk, keQKj, h = 1,...,7, (74) 
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QX^unhi = 2teQB2knih QW^kniht , jsJ, nepR, isQIk, keQKj, h = 1,...,7, (75) 



QX^icnhi = 2t€QB3knih QV^kniiit JeJ, nsFk, isQIk, keQKj, h = 1,...,7, (76) 



5 SisQu Akih Q*ki = I^neFk ^isik QX*knhi , jeJ, keQKj, h = 1,...,7, (77) 



i:ier*kXJki < Dr^" JeJ,kGKj, (78) 



10 



2leQLkQ*Id < QI>*k""^' ,jsJ,kGQKj, (79) 

X^ki, Y^kih, U*kiht, wUt,Vkihi , QX^kriii, Q'ki, QU'knihi, QWW QVkniht ^0 and 
integer for all j, k, n i, h, and t, and G\x> 0 for all j, r, h, and t. (80) 



15 where f Jht(X, Y, Z, U, W, V, QX, Q, QU, QW, QV) is as defined in equation (62) above. 

G^^htin formulation (65-80) above doesn't have to be integer even though it 
represents the number of agents from group j allocated to serve type r contacts. Hence, 
formulation (65-80) is a Mixed Integer Programming (MDLP) formulation involving both 
20 integer and continuous variables. 

Note that MDLP formulation (65-80) may also be modified to have multiple daily 
shift lengths by modifying f \t to include the set of shift lengths Fk, and appropriate akiht 
coefficients for X-^ki, iel , keKj, in (62). 

25 

In order to prevent the notation becoming more taxing, the MILP model of the 
present invention for optimal workforce scheduling in a skills based scheduling 
environment is presented for a combination of the scheduling environments disclosed 
before for formulations (15-23) and (34-41). Depending on the scheduling environment 
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and its tour and shift requirements, the present invention extends the MILP model by 
merging formulations (1-8), (14-23), (24-31), (34-41) or any extension of these as 
disclosed before in a manner similar to the MILP model presented above. It should be 
clear to those skilled in the art that the specific embodiments and the MILP model 
5 disclosed above can be readily utilized as a basis and extended to other scheduling 
environments without departing from the scope of the invention. 

The present invention also extends the MDLP model to obtain a second 
formulation by defining pseudo tours for the tours requiring consistent daily shift start 

1 0 times as follows. Consider a tour type, say tour k, keKj, available for skills group j Set 
Kj contains the valid tour types for skills group j that require consistent daily shift start 
times. Constraints (68-73), (78), and six terms on the right side of (62) involving 
variables X*id, Y^kih, ZV U*kiht, WW and Vkuu specify the break placement and days off 
scheduling requirements and agent availability for tours in Kj. The present invention 

1 5 extends the MILP model for the second formulation by first considering all valid work 
patterns given in FIG. 3 and identifying the patterns satisfying constraints (71-73) for 
keKj. A distinct pseudo tour is then defined for every daily start time specified in Ik, 
keKj, with the vahd work patterns identified. Let one such pseudo tour be a. Then, a 
has only one daily start time in Ola representing one of the daily start times in Ik, has the 

20 same break placement and days off scheduling requirements with tour k, keKj, and has 
the work patterns in OLa . After these modifications, pseudo tours defined this way are 
treated like the tours not requiring consistent daily start times and included in the tours 
set QKj available for agent skill group j. Since Ola has only one daily start time, this 
satisfies the consistent daily shift start time requirement. Also, for tour a. Fa may have 

25 more than one shift length specified for tour k, keKy Another way the present invention 
extends the MELP model is by merging notation used in constraints (68-70) and (74-76), 
and constraints (78) and (79), and removing the redundant notation and constraints as 
follows. 
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Minimize SjeJ SieQU SkeQKj C^U Q^kl + ^jel SkeQKj SreFk SisQIk c'kni QX'kriu (81) 

Subject to 

SjcNir ei^G^'hc > b'ht" ,reR,t8Th, h= ],..., 7, • (82) 

fVQX,Q.QU,QW,QV)- IreNj Gi'h. = 0 , jeJ, teTh, h= l....,7, (83) 

QX'uM = ZteQBiknih QLFkmht , jsJ, nsFlc, ieQIk, keQKj, h = 1,...,7, (84) 

QXW = 2t6QB2iauh QW^knihi , jsJ, neFk, ieQIk, ksQKj, h = 1,...,7, (85) 

QV^knih, , j8J, nsFk, isQIk, keQKj, h = I,... J, (86) 



Siequ Akih Q'ld = 2.KFk SfeDc QXW , jsJ, keQKj, h = 1,...,7, (87) 

SisQUcQ'ki < QD'ic"^ ,j8J,kEQKj, (88) 

QX^kiihi, Q'ld, QU'kniiu, QW^knihi, and QVkniht >_0 and integer for all j, k, n, i, h, t, and 
G'ht>Oforallj,r, h,andt. (89) 

where f \,(X, U, W, V, QX, Y, Z, Q, QU, QW, QV) is now expressed as given below. 

f*h.( QX, Q, QU, QW, Q V) = ZkeQKj I^nsFk Sj^Dc afariht QX^knhi 

- SkeQKj SncFk^iEQTlkiilit QU'kiiiht 

(QW-'k„i,U + QW^knih(t-I)) 
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- SkeQKj 2]nsFk 2iEQT3knht QVkniht 

j8J,tsTh,h=l,...,7, (90) 

Extension of formulations (65-80) and (81-89) to the cases discussed in relation to 
formulations (1-8) and (34-41) is straightforward. It should be clear to those skilled in 
the art that the information disclosed in relation to formulations (l-8_ and (34-41) can be 
readily utilized as a basis to extend formulations (65-80) and (81-89) of the present 
invention to the cases involving (i) minimizing scheduled/paid time or maximizing agent 
preference, (ii) schedule generation when there aren't enough agents (i.e. there are agent 
shortages in some periods for some contact types), (iii) tours with longer and/or more/less 
breaks, (iv) fixed days-ofiF, and (v) days with a closure time earlier than the other days. 
Accordingly, it is to be understood that the embodiments of the invention herein 
described are merely illustrative of the application of the principles of the invention. 
"Referencehereintondetails of theillusfrate^'emb'odfments are not1ntende"d to^rimit^he 
scope of the claims, which themselves recite those features regarded as essential to the 
invention. 

Description of the Method for finding an Optimal Solution of the MBLP 
model and an Optimal Schedule in a Skills Based 
Scheduling Environment 

Formulations (65-80), and (81-89), and their extensions disclosed before can be 
solved to obtain a global optimal solution in a number of ways. Such a global optimal 
solution satisfies the necessary and sufficient conditions for optimality of the MILP 
model of the invention. The solution method of the invention will be discussed in the 
context of formulations (65-80) and (8 1-90). However, the solution method of the 
present invention can also be used to solve optimally formulation (1-8) and its extensions, 
and formulation (24-31) and its extensions as well as formulation (47-59) disclosed 
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before. Therefore, the scope of the solution method of the invention also includes such 
MBLP models. 

The branch and cut (B&C) algorithm is a well known technique for optimally 
solving MILP formulations such as (65-80), and (81-89) (Wolsey, 1998). A computer 
program implementing the B&C algorithm can be developed using any programming 
language or can be obtained from a number of software companies including Ilog, Inc., 
and Dash Optimization. These providers are also offering their B&C algorithm in the 
form of a callable library that can be called from an integrated MILP model generator and 
optimizer module. In either case, the B&C algorithm forms a part of the solution method 
of the present invention. 

The MILP formulation generator ((3) in FIG. 1) of the invention is implemented 
in a computer program. Given the agent and skills requirements for each planning period 
tsQTh, h= 1,. . .,7, and center information regarding agent groups, tours and shift 
parameters, the JVULP formulation generation program develops a formulation (step (11)) 
of the scheduling problem using the MELP formulations and their extensions (1-89) of the 
invention disclosed before. To generate a formulation, the MBLP formulation generation 
program calculates the coefficients of various shift, break, days-off, work pattern, 
overage and underage, agent allocation variables in the respective equations and 
inequalities of the MELP model of the invention for the scheduling environment 
considered. This information is used by the MILP formulation generator program to 
generate data arrays containing the coefficients of the decision variables in every 
constraint defined for the particular scheduling environment to meet its requirements. 
These arrays are then passed to the optimizer module (step (12)). It is to be understood 
that a number of versions of the MELP formulation generation method can be developed 
for the formulations, and their extensions disclosed above without departing from the 
scope of the invention. 

The optimizer module is a computer program working with the MDLP formulation 
generator program, and uses a B&C algorithm supplemented by the RA algorithm of the 
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invention. The purpose of the RA algorithm is to search for integer feasible solutions for 
the MTT.P model generated together with the B&C algorithm. If a sub-problem is 
infeasible (i.e. the sub-problem has no integer or non-integer solution), or has an integer 
feasible solution (that is, a solution satisfying all constraints of the model including the 
5 integrality constraints for all decision variables), or has a solution that is inferior to the 
best available lower bound for the optimal integer solution of the MELP model, the B&C 
algorithm doesn't call the RA algorithm since it can't develop a new or different integer 
solution. Otherwise, the B&C algorithm calls the RA algorithm after obtaining a non- 
integer but feasible solution to a node. At each such node (or after every any number of 
1 0 such nodes), using the non-integer solution found by the B&C algorithm to that node, the 
RA algorithm makes an attempt to find an integer feasible solution to the MELP model of 
the invention while trying to achieve the best objective value. 

The B&C algorithm first solves the continuous version (i.e. Linear Programming 

±5 (LP) relaxation) of an M I LP model ( Wolse y, 1998 ) (ste p (17)). This node i s called node 

zero (subsequent nodes are numbered sequentially). If the continuous solution thus found 
is integer feasible and optimal (step (20)), then the optimizer terminates with an optimal 
solution to the MELP model (step (21)). If this solution is not integer feasible (i.e. if some 
integer restricted variables have non-integer values), then the B&C algorithm starts 

20 creating nodes (sub-problems) by adding constraints on certain decision variables to 
eliminate non-integer values. In this case, the objective value of the solution of the LP 
relaxation is a lower bound on the lowest objective value that can ever be achieved for 
the scheduling environment formulated in the MJLP model. That is, no solution to the 
MELP model can have an objective value better than the lower bound thus obtained. 

25 After solving the LP relaxation at node 0, if the solution is a feasible solution to the LP 
relaxation but violates some integrality constraints (step (22)), the B&C algorithm calls 
the RA algorithm. If the RA algorithm locates an integer feasible solution, its objective 
value is saved as the best heuristic objective value, and the integer solution is returned to 
the B&C algorithmXstep (lOOA)) as the best-known integer solution (if no integer 

30 solution is available beforehand). In the subsequent nodes, the RA algorithm of the 
invention proceeds as outlined below. 
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Step 23: Obtain the objective value for the best integer solution knov^ from the B&C 
algorithm. If this value is better than the best heuristic objective value in the RA 
algorithm, update the best heuristic objective value. 

5 

Step 24: Obtain the values of the shift, breaks, days-ofF, work pattern, and skill 

allocation variables (X, U, W, V, Y, QX, Q, QU, QW, Q V ) in the non-integer 
solution to the node solved from the B&C algorithm. 

10 Step 25: Round the fractional values of all U, W, V, Y, QX, QU, QW, QV, and G 
variables down (i.e. truncate the fractional values). As for X and Q variables, 
round these variables down by their fractional parts are less than or equal to 0.50. 
Keep a list of the variables with a fractional part greater than 0.50. Find the total 
number of agents assigned to each tour type. Round the variables that had a 

.15 fractional part greater than 0.50 u p, one at a time, if the new value is not violating 

the maximum agent constraints (e.g. constraints (78-79)), If the maximum agent 
limit is reached for a tour group, round the remaining variables down even if their 
fractional parts are greater than 0.50. Compute the agent shortages and excesses 
in each period for all contact types using the rounded values of decision variables. 

20 

Step 26: Compare the days-ofF variables Y and Z v^ath shift variables X for tours with 

consistent daily start time and shift length, and daily shift variables QX with work 
pattern variables Q. If additional days off or daily shifts are needed to satisfy tour 
requirements, go to step (27). Otherwise, go to (28). 

25 

Step 27: Schedule more days oflF (i.e. increase the values of appropriate days-oflF 

variables) for the tours with consistent start times on the days with the lowest 
maximum per period, or total shortages during a tour's shift span to satisfy days 
off requirements for every agent currently scheduled. Similarly, schedule daily 
30 shifts to have sufficient number of daily shifts scheduled for the tours not 

requiring consistent daily shift start times on the days and shift start times with the 
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highest maximum per period shortage, or total shortage covered by that shift. 
Thus, the tour requirements for days off, and daily shifts (e.g. constraints (71-73) 
and (77)) will be satisfied by these variables. 

5 Step 28: Compare the break variables, U, W, and V with the shift variables X and days- 
off variables Y, and QX, QU, and QW with daily shift variables QX. If the 
required number of breaks is scheduled for every tour, go to step (30). Otherwise, 
go to step (29). 

10 Step 29: If more breaks are scheduled on a day for a tour, unschedule breaks (i.e. 

decrease the value of appropriate break variables) to have the correct number of 
breaks scheduled. If fewer than the required number of breaks are scheduled, 
schedule additional breaks in the periods (contagious if a break is more than one 
period long within the break window of the respective break) that has the lowest 

15 maximum per period shortage, or total shortage (or if no period has shortage, then 

the highest maximum per period excess, or total excess) to satisfy the break 
balance requirements (e.g. constraints (68-70) and (74-76). 

20 Step 30: Compute the agent shortages and excesses using the agent allocation variables 
G for each contact type, planning period, and day. 

Step 31: Compute the scheduled agent availability in each period for all agent groups 

(i.e. f Jhi(X, U, W, V, QX, Y, Q, QU, QW, QV) in equation (67)), Compare the 
number of agents available in each agent group j in a planning period, say teTh on 
day h, with the values of the agent allocation variables G associated with this 
agent group in planning period t. If the number of agents available to work for 
agent group j is not equal to the sum of the current values of the agent allocation 
variables (from step 25) for all contact types in Nj in period t on day h, adjust 
values of associated allocation variables to meet the requirements by contact types 
in Nj. Prioritize contact types in each period by (i) maximum shortage, (ii) 



25 



30 
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whether this agent group is the only agent group (or the number of agent groups 
and tour groups ) in M\ the set of agent groups that can serve a contact type, (iii) 
lowest cost per period per agent to meet shortage for a contact type. This 
allocation process can be carried out in a number of ways. These allocation 
5 ' criteria are implemented in a nested fashion so that if contact types have the same 

merits as far as criteria (i) is concerned, than the tie is broken by applying criteria 
(ii), and so on. This process results in a schedule that allocates all agents 
scheduled. Update period shortages and excesses. 

1 0 Step 32: The schedule constructed using the non-integer solution found by the B&C 

algorithm for a node satisfies all constraints (e.g. formulation (65-80)), but may 
not satisfy the agent and skill requirements (e.g. constraint (66)). If all agent and 
skill requirements are satisfied in all periods to be scheduled, then go to step 39. 
Otherwise, go to step 33. 

JS ^ 

Step 33: If all available agents are scheduled, go to step 34. Otherwise, go to step 37. 

Step 34: If a feasible schedule satisfying all agent requirements is found either by the 
B&C algorithm, or in the earlier calls to the RA algorithm, go to step 35. 
20 Otherwise, go to step 36. 

Step 35: Discard the current infeasible solution (that is, an integer solution with agent 

shortages). Stop the RA algorithm, and return control back to the B&C algorithm 
(step (17)). 

25 

Step 36: If an integer feasible solution to the MELP model hasn't been found by the B&C 
and RA Algorithms at earlier tries, compare the new infeasible solution with the 
. infeasible solution recorded in earlier calls to the RA algorithm, and update the 
best infeasible solution if the new one has fewer period shortages. If this is node 
30 zero, store the infeasible solution if there are no known solution. Stop the RA 

algorithm, and return back to the B&C algorithm (step (17)). 
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Step 37: Given the current agent allocations and requirements, find the agent requirement 
constraints that are not satisfied (i.e. have shortages) for all contact types and 
periods. Considering all agent groups je J, and tour types keKj and ksQKj, find a 
tour schedule that will reduce the highest number of period shortages per dollar, 
and schedule the tour with an agent. To specify the tour schedule for the new 
agent, find a daily shift start time, days-ofF schedule, and breaks for work days if 
the tour type requires consistent daily start time, or work pattern daily shift start 
times, shift lengths, and break times for work days if the tour type doesn't require 
consistent daily start time and shift length as follows: 

37a: For each agent skill group j with available agents (i.e. the number of agents 
already scheduled in the rounded solution found in step (25) is less than 
the maximum number of agents available for an agent skill group), 

consider tour t y pe k that is available for this a g ent group. If the t our 

requires consistent daily start time and shift length, that is ksKj , go to 
37b. Otherwise, keQKj and go to 37c. 

37b: Consider a daily start times in Pk for tour type k, and all contact types in R. 
Determine the period with the highest shortage level covered by this tour 
(referred to as "peak shortage coverage" in the remainder of this 
document). In order to determine the highest peak shortage coverage by a 
new agent, allocation of an agent's time should be decided based on the 
shortages by contact types to maximize this measure. Given the allocation 
of agent's time, determine also the total number of agent shortage periods 
(referred to as "total shortage coverage" in the remainder of this 
document) covered by this tour. All subsequent references to shortages 
and excesses refer to this allocation of the new agent's time. Consider the 
days-oflF schedule for this agent. Schedule the required number of days off^ 
on the days with the highest total excess during the tour schedule (or if all 
days have shortages, then the days with the lowest total shortage). Days- 
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oflF schedule should adhere to minimum number of consecutive days oflF as 
well as fixed days off requirements. Adjust the total shortages covered by 
this agent. Go to37d. 

37c: Consider a work pattern in QL-'k for tour type k, and all contact types in R. 
For the selected work pattern, determine the daily start time for each work 
day that will result in the highest peak shortage coverage by this tour. In 
order to determine the maximum shortage coverage by a new agent, 
allocation of an agent's time should be decided based on the shortages by 
contact types to maximize this measure. Given the allocation of agent's 
time, also determine the total number of shortage periods covered by this 
tour. All subsequent references to shortages and excesses refer to these 
allocations. Goto37d. 

37d: Consider the work days for the tour identified in 37b or 37c. Schedule the 
relief and lunch breaks for work days into the periods within the specified 
break windows v^th the highest total excess (or, if all periods within a 
break window has shortages, then, then the periods with the lowest total 
shortage). Adjust the total shortage covered by this agent. 

37e: If tour k requires consistent daily shift start times and lengths, compare the 

peak shortage covered by tours starting at different times in Pr and find the 
start time i* with the highest peak shortage covered. If there are ties 
between two or more start times in I\, select the one with the highest 
shortage coverage per dollar 

(total shortage coverage) / CVi 

If tour k does not require consistent daily shift start times and lengths, 
compare the peak shortage covered by tours with different work patterns 
in QL*k and find the work pattern 1* with the highest peak shortage 
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covered. If there are ties between two or more work patterns in QL^k 
select the one with the highest shortage coverage per dollar 

(total shortage coverage) 
(dki H- Zh (number of hours worked on day h) * diaii ) 

37f: If tour k for agent group j is the first tour evaluated, then store the current 
skills group, tour schedule, and agent allocation as the best ones found so 
far, and go to step 37a. Otherwise, if the peak shortage covered by tour k 
is higher than the highest peak shortage covered by any tour and skills 
group considered so far, update the best skills group, tour schedule, and 
agent allocation with the current ones. If there is a tie between the current 
'tour schedule and skills group, and the best skills group and tour schedule 
fdund"s'o"faT,~sele"ct'theio'n^^ 

( total shortage coverage) , for i when keKj, 

(total shortage coverage) 
(Cid + 2h (number of hours worked on day h) * dioii ) 

,for 1 when keQKj, 

AJso update the peak shortage covered, and the total shortage coverage for 
the current tour schedule. Go to 37g if all skills groups and all tours are 
evaluated. Otherwise, go to 37a. 

37g: Add the best tour schedule found in steps 37a-37f to the schedule developed 
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in steps 16-37. Update the values of shift, days-off, break, and agent 
allocation X, Y, U, W, V, QX, Q, QU, QW, QV, and G variables. Update 
shortages and excesses using the expanded schedule with the new agent. 

5 Step 38: If there are shortages in any period and contact type, go to step 33. Otherwise, 
a feasible schedule is found. Go to step 39. 

Step 39: Consider the agents scheduled in the expanded schedule. Start with the most 
expensive skills group and consider the tours scheduled for this group. If there is 
10 a tour that can be removed without creating any shortages in any period and 

contact types, the current solution has redundant tours. Consider all skills groups 
in this fashion to determine redundant agents. Go to step 40 to remove the 
redundant tours. Otherwise, go to step 41. 



1 5 Step 40: Remove all redundant tours together with the shift, break, and days-ofF 
schedules associated with them. 



Step 41: Compare the objective value for the current schedule with that of the best 
feasible solution to the MELP model found so far (if no feasible schedule is 
20 known, then set the cost for best found schedule to a large value such as lO"**^^) by 

both the B&C algorithm and the RA algorithm. Go to step 42 if the current 
solution has a better objective value than the best known feasible solution. 
Otherwise go to step 44. 

25 Step 42: Update the best known feasible schedule with the current solution. Update the 
values received from the B&C for the variables. 
Step 43: Stop the RA algorithm, and return the new solution to the B&C algorithm as the 
current best integer solution (step (16)). 

30 Step 44: Discard the current solution. Stop the RA algorithm, and return to the B&C 
algorithm as the current best integer solution (step (17)). 
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The RA algorithm can also be supplemented with a two-way interchange, or 
general n-way interchange step so that, when an integer feasible solution is found, the RA 
algorithm can make an attempt to improve it by replacing one or more scheduled tours 
with other tours that are not currently scheduled. This is a well-known enhancement in 
the field of optimization that is widely applied to a number of other optimization 
problems. Also, if chosen, the solution method of the invention may also be terminated 
prior to finding an optimal solution if, for example, (i) the objective function for the best 
known integer solution diflFers less than a pre-specified percentage fi'om the best lower 
bound found by the B&C algorithm, (ii) a pre-specified number of sub-problems (e.g. one 
or more sub-problems) are solved by the B&C and RA algorithms, or (iii) a pre-specified 
amount of time elapses after locating the first integer feasible solution to the MDLP. A 
solution that is either optimal or satisfying these conditions is referred to as a "terminal" 
solution. The RA algorithm may also be modified to search for an integer feasible 
solution to a sub-problem using any of the known heuristics (e.g. Tabu Search, Simulated 
Annealing, Genetic Search), or to formulate and add Gomory's cuts or other well know 
cuts (Wolsey, 1998) to sub-problems, change branching direction or priorities in the B&C 
algorithm. The optimization process outlined here doesn't exclude these extensions, 

20 Once a terminal solution is found, detailed tour schedules for agents can be 

developed in a number of ways. All such schedules will have the same objective value 
and, therefore, equally eflfective. First, work patterns and daily shift start times are 
specified. Consider a tour requiring consistent daily start times. The days-ofF scheduled 
in the solution for this tour are first assigned to individual tours. For example, if, in the 

25 terminal solution, tour type 1 has 2 agents fi'om skill group 3 to start daily at 7:00am (that 
is, the variable X'^ktioo) = 2 in formulation (65-80)), initially these 2 agents are scheduled - 
for work on all 7 days in a week. However, since day-ofF constraints (e.g. constraints 
(71-73)) ensures sufficient numbers of consecutive and non-consecutive days off for 
every agent, each of these two agents will get the required number of days-ofF. Suppose, 

30 for instance, each agent assigned to this tour receives only two consecutive days ofF . 
Then, in the solution, Eh Y*'i(7:oo)h = 2. Consequently, each agent will be off on one of 
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these two consecutive days oflF. Thus, the tour schedules for these two agents now have 
work and non-work days specified. Similarly, if a tour does not require consistent daily 
star time and shift length, daily shifts scheduled (including lengths and start times) in the 
terminal solution by the values of the decision variables are assigned to work patterns 
5 specified by the values of the work pattern variables in the terminal solution. Work 
pattern constraints (e.g. constraint (77)) ensures that a daily shift is scheduled for every 
agent that will be at work due to the work patterns scheduled in the solution. These daily 
shifts are assigned to work patterns by following the variable definitions to satisfy tour 
and shift niles for each tour type. Thus, the weekly tours are stated with off days, work 

10 days with daily start times and shift lengths for every agent scheduled in the solution. 

Break balance constraints (e.g. constraints (68-70) and (74-76)) ensure that one break of 
each type is scheduled for every agent scheduled to work on a given day. These breaks 
are assigned to respective tours that have a given day as a work day while making sure 
that agent group, tour type, the daily start time, shift length and break window 

1 5 associations specified in the variable definitions are followed correctly to satisfy break 
scheduUng requirements. Detailed weekly schedules are then assigned to specific agents 
based on their preferences, starting with the agent having the highest preference priority 
(e.g. seniority) to complete the rostering process. 

20 Accordingly, it is to be understood that the embodiments of the invention herein 

described are merely illustrative of the application of the principles of the invention. It 
will be clear to one skilled in the art that the above embodiment may be altered or 
supplemented in many ways without departing fi*om the scope of the invention. 
Reference herein to details of the illustrated embodiments are not intended to limit the 

25 scope of claims, which themselves recite those features regarded as essential to the 
invention. 
What is claimed is: 

1. A computer implemented method for developing an optimized workforce 
schedule for a plurality of agents each with a combination of defined skills to serve one 
30 or more contact types, a plurality of tour and shift scheduling rules located at 
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one or more contact centers each with its own operating hours and time zones, 
comprising the steps of: 

(a) acquiring net staffing level requirements for each contact type, and for each 
period to be scheduled by any means; 
5 (b) acquiring tour, shift, days-ofF, and break scheduling rules, agent skills groups, 

agent availability, and objective criterion to be optimized and its parameters 
by any means; 

(c) formulating the constraints and objective function of a Mixed Integer 
Programming (MILP) model to describe the tour and shift rules, staffing level 

] 0 requirements for a plurality of contact types and for each period to be 

scheduled, agent skills, agent availability, and agent costs; 

(d) solving the LP relaxation of the MILP model formulated in the B&C 
algorithm, and stopping the B&C algorithm if an optimal solution to the MBLP 
model is found; 

15 (e ) callin g the RA algorithm comprising of the following steps_if a solution to a 

node violating some integrality constraints is found by the B&C algorithm: 

i. Obtaining the values of the decision variables in the solution found to 
the MILP model; 

ii. Rounding the fi-actionaJ values of decision variables U, W, V, QX, QU, 
20 QW, QV, and G down, and weekly tour variables X, and work pattern 

variables down if their fractional part is less than or equal to 0.50, and 
up if greater than 0.50, provided the agent availability constraints on the 
maximum number of agents available are not violated; 

iii. Scheduling additional days off^by increasing the values of Y and Z if 
25 additional days off needed to satisfy tour scheduling requirements, and 

additional daily shifts to have a shift scheduled for every agent who will 
be working on a given day based of the work patterns scheduled; 

iv. Scheduling additional daily breaks if the number of breaks of each type 
is not sufficient to satisfy break scheduling rules for a tour, and 

30 unscheduling breaks if more breaks are scheduled due to rounding; 
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V. Computing agent shortages and excesses for each contact type and 
planning period; 

vi. Computing scheduled agent availability for all agent groups, comparing 
them with the agent allocations to dififerent contact types by the rounded 
values of allocation variables G, and adjusting the values of the 
associated allocation variables to make agent allocations equal to 
scheduled agent availability; 

vii. Checking the solution constructed in steps (a) through (f) to find out if 
all agent requirements are met in every period and , when all 
requirements are met, eliminating all redundant agent tours and stopping 
with the integer feasible solution found; 

viii. Stopping with a solution with agent shortages if agent requirements for 
some contact types are not met in all periods and all available agents 
who can serve these contact types are scheduled; 

ix. Continuing to step (i) if there are agents available to schedule and there 

shortages in some periods for some contact types; 
X. Finding all periods in which some contact types have shortages, finding 
an agent with skills needed together with a complete tour schedule with 
work and non-work days, daily shift: start times and shift lengths, daily 
20 break times to reduce the agent shortages, and adding to the solution by 

updating the values of the decision variables to include the newly found 
tour schedule; 

xi. Repeating steps (viii), (ix), and (x) until an integer feasible solution 
satisfying all requirements is found or the algorithm stops in step (viii); 
25 xii. Examining the tours scheduled in the integer feasible solution found and 

eliminating redundant tours 

(f) applying the RA algorithm to the solution found for the current node by the 
B&C algorithm, and passing it to the B&C algorithm if a solution better than 
30 the best integer solution known by the RA.algorithm is found; 
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2. The method of claim 1 fUrther comprising of repeating steps (e) and (f) for • 
every node the B&C algorithm solves and finds a solution violating only some integrality 
constraints, until a terminal solution to the MILP model is reached. 

3. The method of claim 1 further comprising of processing the terminal solution 
found to assign daily shifts with start times and shift lengths to work patterns, and days- 
ofF scheduled to tours, and daily breaks to specific shifts to develop detailed weekly agent 
schedules; 

4. The method of claim 1 fiarther comprising of assigning agents to detailed 

weekly tours; 

5. The method of claim 1, in which a terminal solution to the MILP model 
formulated is found when the objective function value for an integer feasible solution 
differs no more than a pre-specified percentage from the best lower bound found in the 
B&C algorithm for the MILP model. 

6 The method of claim 1, in which a terminal solution to the MILP model 
formulated is found when an integer feasible solution is found and a pre-specified period 
of time is passed. 

7. The method of claim 1, in which a terminal solution to the MILP model 
formulated is found when an integer feasible solution is found and a pre-specified 
number of nodes are solved in the B&C algorithm. 

8. A computer program in a computer readable medium to be run on a computer 
for developing optimal schedules for a plurality of agents with a combination of defined 
skills, a plurality of contact types requiring defined agent skills, a plurality of tour and 
shift scheduling rules, and a plurality of contact centers, comprising of 

means for acquiring agent and skill requirements for each period to be scheduled; 

means for acquiring, for each contact center to be scheduled, operating hours, 

agent skill groups, and scheduling paradigm for contact centers to specify 
whether each center is to be scheduled independently or a plurality of 
contact centers as a virtual contact center with combined; 

means for acquiring tour and shift parameters and scheduling rules, and their 
availability to agent skills and work groups at different contact centers; 

means for acquiring the objective criterion to be optimized and objective function 
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parameters to quantify cost or benefit; 
means for generating an MILP model by defining decision variables, parameters 

and sets, and then formulating the objective function, agent and skill 

requirement constraints, breaks constraints, days-off and work pattern 
5 constraints, agent availability constraints, other schedule related 

constraints, and non-negativity and integrality constraints; 

means for solving the MILP model using an optimization algorithm; 
means for generating a detailed tour schedule using a terminal solution to the 

MILP model found by assigning daily shifts to work patterns scheduled, 
10 - days-off to weekly tours, and breaks to daily shifts, 

means for assigning agents to specific tours. 
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